We claim: 



1. A method for determining areas of structures from an image, the 
method comprising: 

acquiring a digital image; 
. enhancing the digital image using non-linear filters; 

segmenting the digital image using a spatial gradient algorithm, a 
initial point positioning algorithm, and . an , optimal path algorithm that 
combines the spatial gradient and the initial point positioning algorithms to 
produce homogeneous regions; and 

. determining the areas of the hpmog 

2. The method of Claim 1, wherein the non-linear filters include a heat 
filter and a shock filter, the heat filter applied to the digital image followed by 
application of the shock filter to the digital image. 

3. The method of Claim 2, wherein the heat filter is defined by a partial 
differential equation of an inputted image pixel intensity ,u expressed in an equation 
definedas : , 



du __d u d u 



where u is the image being processed. 



dx 



2 . 



is the second partial derivative of u 



along the x-axis, and — r- is the second partial derivative of u along the y-axis. 
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; 4. The method of Claim 2, wherein the shock filter is a partial differential 

' . ■ ' " - ' ' ■ ' ■ ' / ■ 

equation of an imputed image pixel intensity u expressed in an equation defined as 



du . .8 u. 
— = -sign(—) 

ot ax 



du 



dx 



where u is the image being processed, — ^ is the second partial derivative of u, 

dx . 



du 



dx 



is the absolute value of the first derivative of u along the x-axis. 



5. The method of Claim 1, wherein the optimum path algorithm includes 

■ • ■ • ■ ■ ■ ■ ■ ' - ' " " . ' ' ■ '. ' " ■ " ' ' . ■ 

a total cost function C(p) weighted between an edge distance cost function, a path 

direction cost function, and a previous contour distance cost function from the 



10 



equation 



; . C{p) = aC^ (p) + fiC, (p) ^yCXp) . : 



where Cg(/?)is the edge distance cost function, C^(p)is thp path 



direction cost function, Q (77) is the previous contour distance cost 
function, and a^fiyandy are numerical values. ^ , 



15 



6. The method of Claim 5, wherein the edge distance function C^{p) is 



defined from the equation 



P2 



1 



where /j is the pixel intensity, ||V/, is the gradient magnitude the pixel 



at location i, and is a constant value. 
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7. The method of Claim 5, wherein the path direction cost function 
C^ (/?) is defined from the equation 

■ V Q(/o^|;(/;-^r"), \'./- . : 

. ... - t=p\ . ' . 

where is the image intensity of the inside and dark pixel and I^"^ is 
the image intensities the outside and bright pixel along a path bisecting 
the and pixels and connecting pixels points / and U adjacent 

to each side/of the pixel. 

8. The method of Claim 5, wherein the previous contour function C^{p) 
is defined by -the equation , . 



where P is the previous contour and Df is the distance of point i from 
the closest point on the previous contour. 



. 9. A method for determining volumes of structures from a set of images, 

'' ' ■ ' ' ■ ' ' ' -. . " . - ' ' " . ■ ' ■ . ' 

the method comprising: 
15 acquiring at least two digital images; 

enhancing (each digital image using non-linear filters; 

-"■*■ - ^ \ I-"'" • ■** ■ 

segmenting the digital image using a spatial gradient algorithm, a 
initial point positioning algorithm, and an optimal path algorithm that 
combines the spatial gradient and the initial point positioning algorithmto 
20 produce homogeneous regions; 
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assembling the digital images into an array; 

determining the areas of and the volumes between the homogeneous 
regions in the array. 

10. The method of claim 9, wherein the array includes a rotational 
assembly, a wedge assembly, and a translation^ assembly. ^ 

11. The method of Claim 9, wherein the non-linear filters include a heat 

filter and a shock filter, the heat filter applied to the digital image followed by 

• ' . ' ■ - ' ■ ' " - ' " ■ • ■ ' ■ ■ , ' .■ " • , 

application of the shock filter to the digital image. 

. 12. The method of Claim 11^ wherein the heat filter is defined by a partial 



10 



differential equation of an inputted image pixel intensity u expressed in an equation 



defined as 



where u is the image being processed. 



is the second partial derivative of u 



15 



along the x-axis, and -^-y is the second partial derivative of u along the y-axis. 



13. The method of Claim 11, wherein the shock filter is a partial 



1 ■ t 



differential equation of an imputed image pixel intensity u expressed in an equation 
defined as 



du 
'dt 



ax 



du 



dx 
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where u is the image being processed, — i- is the second partial derivative of u, 



du 



dx 



is the absolute value of the first derivative of u along the x-axis. 



14. The method of Claim 9, wherein the optimurn path algorithm includes 

' • ' ■ . -. ■ ' . . ■ ' ■ ^ ■ -■ . ■ . ■ 

a total cost function C{p) weighted between an edge distance cost function, a path 



5 



direction cost function, and a previous contour distance cost function from the 



equation 



where Q (p) is the edge distance cost function, . (/?) is the path 
direction cost function,Q.(/7)is the previous contour distance cost 



10 



function, and a,j3,and/ are numerical values. 



15. The method of Claim 14, wherein the edge distance function C^ip) is 
defined from the equation . ^ 



) 



where // is the pixel intensity, |V/, is the gradient magnitude the pixel 



15 



at location i, and // is a constant value. 



16. The method of Claim 14, wherein the path direction cost function 
C^(p) is defined from the equation 



p2 



c,{p)=^ir-ir)^ 
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where //" is the image intensity of the inside and dark pixel and //^^ is 
the image intensities the outside and bright pixel along a path bisecting 
the and //'"^ pixds and coimecting pixels points / and adja^^ 

; to each side of the //'' pixel. 

\ 17. The method of Claim .14, wherein the previous contour function 
C (/?) is defined by the equation , 



P2 



' • /-pi . ■ 



where P is the preyious contour and Z)f is the distance of point i from 



the closest point on the previous contour. 



18. A method to determine volume of a structure in digital images 

> ^ • . . t . " - • 

acquired from electromagnetic and non-electromagnetic sources, the method 



comprismg: 



positioning a transceiver exterior to a patient such that at least a 



portion of the structure is within a field of view of the transceiver, the 
transceiver configured to send electromagnetic radiation and to receive echoes 
of the electromagnetic radiation; 

sending the radiation from the transceiver to the structure; ^ 

■ » . 1 ' ' • " - 1 * 

receiving echoes of the radiation reflected from the structure to the 
transceiver; 

associating the received echoes to form a plurality of 2D scanplahes so 
that they form an array; 
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enhancing the images of the structure in each plane of the array using 
non-linear filters; and 

determining the structure volume spanning between and through each 
plane in the array: 

19. The method of Claim 18, wherein plurality of 2D scanplanes are 
assembled into a plurality of arrays including a rotational array, a translational 
array, or a wedge array. - 

20. The method of Claim 18, wherein the non-linear filters include a heat 
filter and a shock filter, the heat filter applied to the digital images followed by . 
application of the shock filter to the digital images. 

21. The method of Claim 20, wherein the heat filter is defined by a partial 
differential equation of an inputted image pixel intensity u expressed in an equation 
defiinedas 

„ /Where u U .h. u^bei^ p^c^seA, |^ is «,e second paHM derivative of 



u along the x-axis, and is the second partial derivative of u along the y- 



' axis. . ■ . •- 

' - • . . ' * . I . . - ^ . - 

t ^ \ • " - • . ' . . ' 

22. The method of Clmm 20, wherein the shock filter is a partial 
differential equation of an imjjuted image pixel intensity u expressed in an equation 
defined as 
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dt 



dx 



du 



dx 



where u is the image being processed 



is the second partial derivative of u, 



du 



dx 



is the absolute value of the first derivative of u along the x-axis. 



23. The method of Claim 18, wherein the plurality algorithms includes a 
spatial gradient algorithm, a previous contour algorithm, and , an optimal path 
algorithm that combines the spatial gradient and the previous contour algorithms to ; 
"produce the homogeneous regions. . 
■ 24. The method of Claim 23, wherein the optimum path algorithin 
includes a total cost function C(/7) weighted between an edge distance cost 



10 



function,, a path direction cost function, and a previous contour distance . cost 



function frorn the equation 



C{p) = aC,{p) + pC,{p) + yCXp), : ' . : 

where C^{p) is the edge distance cost function, (/?) is the path 

direction cost function, C_ (p) is the previous contour distance cost 



15 



function, and a,p,andy sltq numerical values. . 



25. The method of Claim 24, wherein the edge distance function C^{p) is 
defined from the equation ■ 



. A?2 1 

Ce(p) = £(— : 
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where // is the pixel intensity, V/. is the gradient magnitude the pixel 



at location i, wad ju isconstant value. 



['.■"■ ''-: 26. The method of Claim 24, wherein the path direction cost function 
Cj(p) is defined from the equation 



p2 



V 



where //" is the image intensity of the inside and dark pixel and is 
the image intensities the outside and bright pixel along a path bisecting 



the /-'" and I^""^ pixels and connecting pixels points / and i,i adjacent 



10 



to each side of the I '" pixel. 



27. The method of Claim 24, wherein the previous contour fimction 
C^{p) is defined by the equation . ■ 



p2 



c.(p)=XA^ 

i=p\: 



where\P is the previous contour and Z).^ is the distance of point i from the 



15 



closest point on the previous contour. ; 
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